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Abstract 

For a deconfining thermal SU(2) Yang-Mills plasma we discuss the role of 
(anti)calorons in introducing nonthermal behavior effectively parametrized in 
terms of Planck's quantum of action h. While nonthermal features in thermo- 
dynamical quantities cancel exactly between the ground-state estimate and 
its free quasiparticle excitations they resurface as small radiative corrections 
at a fixed loop-order. Next, we consider thermal 2— > 2 one-loop scattering 
of those modes that remain massless upon the (anti)caloron induced adjoint 
Higgs mechanism. In a stepwise procedure, starting with stringent analytical 
arguments and concluding with a Monte-Carlo based analysis, we argue that 
the constraints, imposed by the emergent, inert, and adjoint scalar field 0, 
are sufficiently strong to render photon-photon scattering within the Yang- 
Mills gas void of measurable consequences. By the optical theorem, we thus 
conclude that pair creation/annihilation of massive vector modes by/into two 
photons cannot be observed. An adaption of the here-developed methods to 
the analysis of irreducible bubble diagrams may prove the conjecture of hep- 
th/0609033 on the termination of the loop expansion of thermodynamical 
quantities at a finite irreducible order. 



1 Introduction 



Effective, indeterministic behavior, inherent to squared scattering amplitudes in 
quantum Yang-Mills theory, could be associated with the presence of classical field 
configurations of finite action and nontrivial topology in 4D Euclidean spacetime 
[U [21 El H] • For SU(2) Yang-Mills thermodynamics a concrete argument supporting 
this idea was put forward recently E] . Namely, in its deconfining phase an ad- 
joint Higgs mechanism, caused by calorons and anticalorons of topological charge 
modulus unity, which effectively materialize as an inert scalar field upon spatial 
coarse graining [7J, splits the mass spectrum into two modes always propagating on 
their quasiparticle mass-shell (unitary gauge) and a massless mode (photon) whose 
off-shellness is constrained by \<p\ (Coulomb gauge) [8]. Moreover, (anti)caloron me- 
diated 4-vertices are pointlike in the effective theory because four-momentum trans- 
fer through them is restricted by \<f>\ in all Mandelstam variables. The present paper 
intents to investigates the consequences of these constraints for photon-photon scat- 
tering in the effective theory for the deconfining phase. However, before performing 
the according analysis, we would like to elaborate on the argument for and conse- 
quences of (anti)caloron mediation of effective Yang-Mills vertices. To do this, we 
work in units where the c, speed of light in vacuum, is set to unity but Boltzmann's 
constant ks and Planck's quantum of action h are dimensionful. 

Due to Feynman [HI OH] and Schwinger [TT] the fundamental Yang-Mills partition 
function Z is representable in fundamental variables as follows 

<7 — 4. -Y-K S^J I 4- -TT-^I \ f Tl/1 - -r~ dr' d z x Ce^} n , 

Z = tr e k b = 2_,\ n \ tr e fc s \n) = VA e k s J0 , (1) 

n ' .M M (0,x)=A M (/3,x) 

where the functional integration is over gauge-inequivalent periodic gauge-field con- 
figurations, subject to an exponential weight associated with their classical action, 
see also text below Eq.|4]). 

A simplification in performing partition function Z in Q occurs if the quantum 
states \n) are taken to be (gauge-inequivalent) energy eigenstates with eigenvalues 
E m of the quantum Yang-Mills Hamiltonian "H. Formally, the partition function 
then reads 

M(m)e k B . (2) 

m 

Here M(m) denotes the degeneracy of E m . Resorting to a representation of the 
Hilbert space in field variables, % is a gauge-invariant functional of the gauge field 
and its conjugate momentum n^. The latter's quantization conditions (equal- 
time commutators of A^ and n M ) involve h. In a Euclidean spacetime, the operator 

e k B formally evolves a state \n) in r - a variable of dimension inverse temperature 
T - from t' = to r' = /3 = 1/T if r' and time r are related as 

, k B h 
t = —t or r = - — r (6) 
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Now, neither the macroscopic concept time (t) nor the macroscopic concept tem- 
perature (t') ought to depend on the values of h and ks- As a consequence, Eq. ^ 
implies that in the formal limit h — > also ks — > such that 

h 

- — = const . (4) 
k B 

Quantity (n|e fc s \n) thus represents the formal persistence amplitude for state 
|n) under (Euclidean) time evolution from r = to r = ^Ap. In field space, 
the summation of this amplitude over \n) can be cast into a functional integral 
over gauge-inequivalent, periodic field configurations subject to measure VA 
and weight e fe s ^° dr d xC eI a ^}^ gee p ar ^ Q f jr;q_ (]_) w here denotes the 

classical Euclidean action density of configuration A^(r, x) after the substitution 
r = t^t'. 

k B 

If the world were purely thermal and Euclidean (no reference to a Wick rotated 
version of time evolution acting on a quantum state in Minkowski spacetime) on 
all subatomic length scales then the concept of action would be alien, and energy 
per temperature in units of ks would be natural to occur in exponential weights 
representing the Yang-Mills energy spectrum. Conversely, if the subatomic world 
were purely Minkowskian and nonthermal then one would work with the action of 
isolated field 'trajectories' in units of h to describe temporal evolution of quantum 
states in terms of the exponentiated Hamiltonian %. Thermal Quantum Yang-Mills 
theory as a whole emerges as an interplay of these two concepts. For example, the 
notion of the thermal ground state invokes the fact that the former concept can be 
related to the latter if effective variables are (nonlocally) represented by fundamen- 
tal ones: The periodicity of these fundamental (classical, Euclidean) configurations 
(calorons and anticalorons) is a consequence of identification (|3]) which invokes time 
evolution of quantum states. Collectively, however, nontrivial periodicity collapses 
to a mere choice of gauge [8]. On the other hand, the occurrence of effective quan- 
tum corrections, expanded in powers of h, appears to be related to the presence of 
(anti) calorons. Notice that the physics of ( ant i) calorons can be interpreted classi- 
cally hinging on adiabatic slowness to render Euclidean signature irrelevant [T2l H~3] . 

For the deconfining phase of SU(2) Yang-Mills thermodynamics we have in 016] 
identified the Euclidean action of relevant calorons or anticalorons [1~4"1 [T5l [T6| [T7] 
with h. The according chain of arguments relies on a Minkowskian interpretation 
of free massless or effectively massive (adjoint Higgs mechanism) one-loop fluctu- 
ations in thermodynamical quantities about the ground-state estimate. Here the 
contribution of these (quasi)particles to the pressure, expressed by a Matsubara 
sum in Euclidean spacetime, is cast into a continuous integration by employing 
Cauchy's integral theorem. Because of Minkowskian on-shellness (inertness of the 
adjoint scalar field <fi [12]) only the Bose weighted part in the massive propagator 
contributes to the loop integration. Thermodynamic consistency (linking of ther- 
modynamical quantities by Legendre transformations motivated by the left part of 
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Eq. (jl|) and a dimensional counting in effective interaction monomials with appro- 
priately normalized field variables [18] determines the effective coupling e in units of 
fir 1 / 2 [HIE]. Finally, one appeals to the observation that (anti)caloron radii p ~ |</>| _1 
dominate the spatial coarse-graining - generating the field cf> [19] - to justify the 
use of e in the formula 5* = for the Euclidean action of the fundamental field 
configurations caloron or anticaloron with topological charge modulus unit}Q As 
an immediate consequence, pointlike effective vertices, introducing perturbations 
into the orderly, thermally weighted Minkowskian propagation of free plane waves, 
occur because these vertices effectively interpolate the incoming and the scattered 
wave through spatially unresolved (anti)calorons subject to a Euclidean time de- 
pendence. Here, the irreconcilability of Euclidean and Minkowskian time evolu- 
tion within a vertex represents the cause of effective indeterminism in the outcome 
of the scattering event: Wick rotation of a caloron's time dependence generates 
a field configuration which is far away from the nearest classical solution to the 
Minkowskian Yang-Mills equations. Effectively parameterized by h, the action of 
isolated (anti)calorons thus introduces disorder to thermodynamically equilibrated 
Minkowskian plane- wave propagation]^] In this sense (anti)calorons behave nonther- 
mally, resolving the following paradox associated with the ground-state estimate: 
While the latter, describing the effects of interacting (anti)calorons in terms of the 
field and a pure-gauge configuration a gs after spatial coarse-graining, obeys an 
equation of state P gs = —p gs it does exhibit a finite heat capacity. In other words, 
the ground state's entropy density s gs , given as 

s gs = ^(P gs + p gs ) , (5) 

vanishes although p gs = 4nTA 3 is temperature dependent, A the Yang- Mills scale 
[8]. At finite temperature this would violate the third law of thermodynamics. 
Therefore, a thermodynamical interpretation of the ground-state physics alone is 
not admissible. That the ground-state estimate does not behave thermodynamically 
by itself is, however, no longer a surprise because of the de-thermalizing effect of 
the ground state's constituent (anti)calorons visible in isolated, effective scattering 
events. The overall Yang-Mills system, however, behaves thermodynamically due 
to the T dependence of the effective coupling e which cancels nonthermal effects 
between the ground state and its quasiparticle excitations [20] . When quasiparticles 
disappear entirely at the critical temperature T c / (lowest attainable temperature 
in the preconfining phase [H]) the thermal ground state would represent the entire 
thermal partition function. This does not make sense. Indeed, the onset of the 
Hagedorn transition at T c > inhibits a thermodynamical description of the Yang-Mills 

^t was shown in [7] that (anti)calorons of higher charge modulus do not contribute to <fi. 

2 This equilibration is by interaction with the ground-state estimate which, in turn, invokes a 
sum over an infinite series of effective forward-scatterings (Dyson series) induced by (anti)calorons: 
the emergence of quasiparticle mass by the adjoint Higgs mechanism. Notice that the Bose-Einstein 
function does not depend on h separately but only on the ratio ^p, see discussion of Eq. (W|) above. 
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system [20J. Except for the limit T \ T c |8], where radiative corrections vanish [8], 
nonthermal effects are ubiquitous also in the deconfining and preconfining phases: 
there are small radiative corrections associated with fixed or resummed loop orders 
[2T| 122] 123] . triggered by isolated (anti)calorons, in the former [21], and there is 
supercooling subject to ground-state tunneling in the latter phase [2D] . 

This paper is organized as follows. In Sec. [2] we review the effective theory to- 
gether with its constraints in unitary- Coulomb gauge. Sec. [3] discusses how these 
constraints apply to exclude certain one-loop diagrams potentially contributing to 
the photon-photon scattering amplitude. Moreover, a technology is developed that 
allows to systematically discriminate certain channel combinations for the two 4- 
vertices in the admissible diagrams contributing to this amplitude. For the remain- 
ing possibilities, where no further analytical assessment can be made, a Monte-Carlo 
simulation of the domain of integration is performed and selected hits are investi- 
gated further in view of local properties of the extremely filamented algebraic vari- 
eties comprising the domain of integration. Finally, we analyze to what extent Bose 
suppression of internal lines influences the result for the amplitude. The last section 
summarizes our results and gives conclusions. 

2 Effective theory, notational conventions, and con- 
straints 

For the benefit of the reader we revisit briefly the effective theory for deconfining 
SU(2) Yang-Mills thermodynamics and discuss its constraints on 4-vertices, see also 
[6j [19] . From now on we work in supernatural units H — ks — c — 1. 

2.1 Effective action 

The effective action for the deconfining phase emerges upon a spatial coarse-graining 
over free, trivial-holonomy (anti)calorons of charge modulus unity and plane- wave 
fluctuations |7J. While the topologically nontrivial sector of field configurations 
yields an inert, adjoint scalar field coarse-graining over interactions within the 
topologically trivial sector is determined by perturbative renormalizability [^} The 
physics of this sector contributes to the effective action for the coarse-grained, topo- 
logically trivial gauge field a M in the same form as the fundamental Yang-Mills action 
[211 [25]. Moreover, the appearance of mixed operators of mass dimension higher than 
four, involving both fields a M and 0, is excluded because momentum transfer to the 

3 Because, a priori these interactions between fundamental modes are characterized by momen- 
tum transfers > \(f>\ 2 loop momenta are far off their mass shell. Thus reference to thermodynamics 
is required in a formal coarse-graining over the topologically trivial sector. Moreover, the sub- 
traction of ultraviolet divergences in perturbation theory is consistent with the fact that calorons 
of radius p ~ |^| sharply dominate this coarse-graining. Thus only a thin shell in momentum 
transfer in fundamental plane- wave scattering is effectively mediated by them [3 [S] . 
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field is impossible [T9"] . Thus the unique (Euclidean, gauge-invariant) density of 
effective action Asfa^] is given as 

£ cff [a M ] = tr - G^G^ + (D^f + — , (6) 



where = d^dy— d u a^— ze[a M , a v ] = £ a denotes the field strength, 1)^0 = <9 M </>— 
ie[a^,0], and e is the effective gauge coupling yet to be determined. £ efl in Eq. ^ 
yields a highly accurate tree-level ground-state estimate and, as easily deduced in 



unitary gauge = 2\(j)\ t s , a tree-level mass m = 2e\(j)\ = 2e^j for propagating 

gauge modes a*' 2 . Because the mass of these modes is generated by a summation 
of a Dyson series of up to infinitely many local interactions with the field <fi 2 the 
off-shellness inducing absorption or emission of propagating modes is forbidden: 
these processes would introduce momentum transfer to the inert field <p. As a 
consequence, massive modes propagate thermally on their mass shell [19]. Massless 
modes (photons), on the other hand, may be off their mass shell at most by \<j)\ 2 [8]. 
Finally, thermodynamical consistency of the system's ground-state estimate and its 
free quasiparticles yields the value e = VEtc almost everywhere. The exception is a 
narrow, logarithmic pole at T c = 13.87^ [8]. 

Considering radiative corrections, the effective 3- and 4-vertices are given as 

r^ c = e(2v)%p A + p B + p c )e abc x [g^ (p B - p A f 

V 7 (PC ~ VbT + 9 ia (Pa - pcf] , (7) 

rffih = ~ ^ 2 (2^5(p A +p B +pc+ Pd) x [6^(sY - g a5 g M ) 

+e ace e bde (g^ g' 5 - g a5 g^) + e ade e bce (g^ 'g* ~ sTg")] , (8) 



compare with Fig.[TJ In the present paper, no explicit reference to propagators is 
made except for the above mentioned fact that massive modes always are on-shell. 
Let us now set up our notational conventions. 

Vertices. Vertices are discussed on different levels. Some arguments require only 
formulas that are valid independently of the sort of attached gauge modes. In such 
cases we just use the symbols p±,, P2 etc. for four-momenta. On a level, where we 
would like to distinguish between massive and massless modes, we use R, S, P for 
the four-momenta of the massive modes and p, q for the four-momenta of photons. 
In applications to the actual scattering process we employ a and b to indicate the 
four-momenta of incoming photons. Outgoing photons are labeled by c and d. Four- 
momenta of internal, massive (loop) modes are denoted by u and v. 

Feynman diagrams. In general, the propagation of massive particles is indicated by 
double lines in Feynman diagrams while massless modes are represented by single 
lines. 
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Figure 1: Yang-Mills vertices, px being a four-momentum, greek letters are Lorentz 
indices, and lower case latin letters denote Lie- Algebra indices. 

Scattering channels. An overall channel in photon-photon scattering is labeled by its 
Mandelstam variable S,T, or U (capital letters). On the other hand, for scattering 
channels associated with a given 4-vertex Mandelstam variables are in lower case 
letters. In this way, one specific configuration can be written in a short-hand nota- 
tion. For example, "Stu" describes overall S-channel scattering with the t-channel 
realized at the first 4-vertex and the u-channel at the second 4-vertex. 

2.2 3- vertex 

At the 3- vertex, given in Eq. ([7]), a photon alway connects to two massive modes 
Let us now check whether the photon can be an on-shell particle of positive 
energy connecting to a 3- vertex Because of four-momentum conservation and the 
on-shellness of the massive modes the following conditions apply: 

I? = S 2 = m 2 , 
p 2 = 0, 

(R + S) 2 =p 2 = 0. (9) 

The energy of the massive particles can be positive or negative, depending on the 
direction of energy flow in the loop of the overall scattering diagram, see Fig. [3} 
Therefore, one has 

4 This is due to the Levi-Civita symbol e a bc in Eq. ^ and expresses the conservation of unbroken 
charge in U(1)gSU(2). 

5 In the present paper we neglect modifications of the free dispersion law po — |p| due to a 
resummation of insertions of the one-loop polarization tensor, see [211 E2] ■ These effects are sizable 
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Figure 2: Two different massive modes (dashed and dotted) and a massless one 
(wavy) at a 3-vertex. 





Figure 3: Box diagram (left) and penguin diagram (right). Both possibilities are 
excluded to contribute to photon-photon scattering. 



IRI 



m 



2m 2 ± 2y/\R\ +m 2 y |S| 2 + m? - 2RS = 
') (ISI 



m 2 ) 



(RS) 



2m 2 RS = 



m 



|R| 2 + |S| 2 + 2|R| |S|cosZ(R,S)) 

+ |R| 2 |S| 2 (l -cos 2 Z(R,S)) = 0. 



(10) 



The minimal value of each summand in the last equation of ( 10 ) is zero: The first 
summand vanishes for cos Z(R, S) = — 1 and |R| = |S|, the second summand is zero 
for cosZ(R, S) = ±1. So the only configuration satisfying Eq. (10) is Z(R, S 



7T 



and equal energy of the photons. Since in a loop integral the thus allowed integration 
over angular variables is over a hypersurface of measure zero we conclude that 3- 
vertices do not contribute to the overall one-loop scattering of photons, compare 
with Fig. [3} 



at low energy and temperatures, and they die off in an exponential way in the former and a 
power-like way in the latter variable. 



2.3 Constraints on 4- vertex 



It is easy to show that the 4-vertex is invariant under a permutation of the legs 
attached to it [26J. This implies that the 4-vertex does not distinguish between 
Mandelstam variables s, t, and u in mediating 2—7-2 scattering. The according con- 



straints in the effective theory, however, do. We have [H] 

S= \{Pl+P2?\ = |(P3 +P4) 2 | < \<P\ 2 , (11) 

t=\(Pl-P3) 2 \ = \(P2-P4) 2 \ < |0| 2 , (12) 

u = \ {p 2 -p 3 ) 2 | = |(pi — P4) 2 | < |0| 2 (13) 



for a scattering amplitude mediated by a 4-vertex with momentum labeling as de- 
fined in Fig. [4] 




Figure 4: Scattering due to a 4-vertex. 

It is also straightforward to demonstrate that, modulo leg permutation, only the 
two diagrams depicted in Fig. [5] contribute. 




Figure 5: Photon-photon to massive-massive and massive-massive to massive- 
massive scattering. These are the cases that need to be distinguished in all scattering 
channels. 
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Specializing Eqs. (11), (12), and (13) to the situation of Fig. IB] and introducing x G 



{p, P} to cover both diagrams, we have 

s = \(R + S) 2 \ < |0| 2 
t = \( X -R) 2 \ < \<f)\ 2 
' (x - S) 2 \ < \<P\ 2 



u 



(14) 
(15) 
(16) 



Let us now show that it is impossible to satisfy all of these constraints simultane- 
ously. Appealing to Eqs. (14) through (16) and using P 2 



R 2 = S* 



m 2 , p 2 = 0, 



one arrives at the following constraints for the case of two photons, two massive 
modes (left) and four massive modes (right): 



two photons, two massive modes: 



four massive modes: 



8= \2R S + 2m 2 -2RS| < 



(17) s = \2R S + 2m 2 -2RS| < |0| 2 , (20) 



2p R + m 2 + 2pR\ < |0| 2 , (18) t = \-2R P + 2m 2 + 2RP\ < 



(21) 



u 



2p 5'o + m 2 + 2pS| < 



(19) 



u 



25 , P + 2m 2 + 2SP| < |0| 2 . (22) 



We have 



\2p R \ > \2p U R 2 - m 2 ) cos (Z (p, R)) | = |2pR| , 



(23) 



|2i? 5 | > \2{^R 2 - m 2 ^S 2 - m 2 ) cos(Z (R,S))| = |2RS| . (24) 
The following argument does not rely on a restriction of the signs of the energies po, 



Ro, and So- From (24) and the fact that m > |0| it follows that the sign of both R 



and Sq in (17) need to be different. Similar conclusions can be drawn for the other 



five cases (18) through (22). Again using variable x G {p, P} we can summarize the 
situation as 



s < 
t < 
u < 



| 2 sgn( J R ) = -sgn(So) 
| 2 sgn(x ) = sgn( J R ) , 
| 2 sgn(xo) = sgn(S'o) . 



Therefore, one arrives at the following statement: 
sgn(xo) = sgn( J R ) = -sgn(So) 



-sgn(xo) . 



(25) 
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As a consequence, the constraints (17) through (19) or (20) through (22) can not be 



satisfied simultaneously. Thus, rather than imposing all constraints simultaneously 
at a given 4-vertex, a coherent superposition of the contributions of nontrivial scat- 
tering channels needs to be considered. As mentioned above, the every channel is 



treated in the same way. For the case, where none of the constraints (18) through 



(22) is trivial, this leads to the following prescription for the implementation of a 



4-vertex in the effective theory [B]: 



[A]abcd 



[4]abcd 



+ r 



afi'yS 
[i]abcd 



+ r 



afS-yS 
[A]abcd 



(26) 



If, on the other hand, say, the t channel is trivial (polarization tensor of the massless 
mode on the one-loop level or figure-eight two-loop contributions to the pressure 
[2T| [22] ) then the vertex acts as 



[4]abcd 



[A]abcd 



+ r 



[i]abcd 



(27) 



This decomposition of a given 4-vertex into its scattering-channel contributions is 
the reason for the large number of 3 3 = 27 possible scattering-channel combinations 
in the three overall channels for one-loop photon-photon scattering, compare with 
Fig.[6j 

d. 



u 



-Ch 





Figure 6: The overall scattering channels for the only admissible diagram in one- 
loop photon-photon scattering. Here a,b,c, and d are the four-momenta of the pho- 
tons, and u and v denote the loop four-momenta of the massive modes. 



2.4 Energy- flow constraints 

Let us now investigate all possible constraints that can occur in one-loop photon- 
photon scattering in view of energy flow. To do this, it is advantageous to work 

6 It is easy to show that the integrand in this case is invariant under k — > — k (k a loop four- 
momentum) which renders the contributions of s and u channels equal. Thus it suffices to compute 
the s channel contribution with weighting unity as it was done in [211 122] . 
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with dimensionless quantities. We will use two different normalizations. A tilde 
on top of a four-momentum component x marks normalizations with respect to the 
modulus of the adjoint scalar field \<p\, a hat denotes normalization with respect to 
the temperature T: 



x 



2e 



x 



-x 



A2 X 

2^f 



A^ 



m 27i 1 2tt 
Here the dimensionless temperature A is defined as 

2nT 

A 



-x . 



A 



(28) 



(29) 



(30) 



where e > \/8n [20]. 

For a 4- vertex connecting two photons (four-momenta p, q) with two massive modes 
(four-momenta R, S) simple combinatorics allows the following possibilities, redun- 
dant by four-momentum conservation and constrained by the maximal resolution 
|0| of the effective theory, for momentum transfer through the vertex: 



Furthermore, the dimensionless mass m = rS- is given as 

m =2e , 



~\ 2 



~\ 2 



R±S 

p±R) 

p±S 

q±R 

q±S 



< 1 

< 1 

< 1 

< 1 

< 1 



\(p±q) 2 \ < 1 



Both photons and massive modes are on shell: 

Ro 



± 



R 



+ m 



> t - 



So 


= ± 


P 2 


= 0, 




= 0. 



+ m 



(31) 

(32) 

(33) 

(34) 

(35) 
(36) 

(37) 

(38) 

(39) 
(40) 
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From Eqs. (30) and (37) through (40) it follows that certain sign combinations of 



Ro and Sq are forbidden. Let us now classify those excluded cases. To do this, we 
consider 



1 > 



_A 2 

R±S 



2m 2 ± 2R S =F 2RS 



2m 2 ± 2 ± 



R 



+ m 2 ± 



+ m 2 =f 2RS 



(41) 



and 



1 > 



p±R 



m 2 ± 2p i?o =F 2pR 



m 2 ± 2p I ± 



R 



m 2 =f 2pR 



(42) 



In the following discussion cases (31) are treated in terms of (41) while cases (32) 
through (35) are covereo^Jby (42). To proceed, notice that 



± 



R 



+ m 2 \ ± 



m? 



r,2 



^RS>m 



(43) 



and 



Po ± 



R 



m 2 



TpR > o 



(44) 



vectors 



R 



Inequality (|43|) is true because of the Cauchy-Schwarz inequality applied to the two 
i \ and ( S , rh \ with the canonical scalar product of IR 2 and the fact 

< 1. Inequality (44) is selfevident. From (43), m 2 > 1 (see (30)), 



that 



cos ZRS 



and (41) it follows that 



sgn ( Ro ) = ±sgn ( S { 



are forbidden for cases (31), respectively. Because of (44), m 2 > 1, and (42) it is 
clear that 

sgn (r ) = ± 



are forbidden for cases (32) and (34), respectively. In addition, the respective cases 



(33) and (35) (R replaced by S) are also excluded. Obviously, no implication for 



the signs of R or Sq arises from ( 36 ) . 



7 We do not distinguish R and S or p and q. For the argument it is only important that the 
former are associated with massive modes and the latter with photons. 
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3 Photon-photon scattering 



Let us now assess the options for one-loop photon-photon scattering. To do this, we 



first apply the results of Sec. 2A to each of the overall scattering channels depicted in 
Fig. |6j This excludes a vast majority of a priori thinkable energy-flow combinations. 
In the overall S-channel one of the remaining six scattering channel combinations 
for the two 4-vertices can be excluded by further analytical arguments. All other 
yet possible combinations are investigated by a Monte Carlo like numerical analysis. 

At a given vertex, we use the following notation to keep track of excluded com- 
binations of the two loop energies Uq and vq. Every cell in a 2 x 2 table accounts for 
one possible combinations of energy flow: 



uo > 0; Vq > 


MO > 0; Vq < 


Mo < 0; Vq > 


Uq < 0; Vq < 



An "X" in a given cell signals that the corresponding combination is forbidden. 



3.1 Energy flow: Overall S channel 



Vertex 2 




Vertex 1 



Figure 7: The overall S-channel. 



The lower vertex in Fig. [7] is denoted by number 1, the upper one by number 2. At 
a given vertex the constraint on a scattering channel can be expressed in a twofold 
way because of total four-momentum conservation across the vertex. For example, 

~\2| 



the constraint at the first vertex 



~\ 2 
a + b 



< 1 can be rewritten as \ (u — v) \ < 1. 



To exclude all sign combinations of the energies in the loop momenta one needs to 

look at either form and combine both statements 

Let us now visualize the constraints at vertex 1 in terms of their tables: 



8 The two 2x2 tables, obtained from each of the forms expressing the vertex constraint, are 
put on top of one another. 
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a + b) 



s-ch: 1 > 



t-ch: 1 > I (a — u) 2 \ 



u-ch: 1 > I (a + v) 2 \ 



(u — v) 2 \ —> 





X 


X 





For vertex 2 we obtain: 



s-ch: 1 > \(u — v) 2 \ 



b + v 



b — u 



~\ 2 
c + d 



X 




X 


X 



X 




X 


X 



t-ch: 1 > 



u-ch: 1 > 



u 



~n2| 
C 



~\ 2 

v + d 





X 


X 





X 




X 


X 



u — d 



\{v + cf \ -)■ 



X 




X 


X 



In the diagram for the overall scattering channel, see Fig. [7j the constraints on the 
two vertices have to be satisfied together in order to contribute: If one combination 
of energy flow and scattering channel combinations is forbidden by one constraint 
then it is not allowed at all. Thus it is suggested to use a 3 x 3 table for visualization 
of excluded scattering-channel and energy-flow combinations within a given overall 
channel. Each of the cells of this 3x3 table, corresponding to a certain combination 
of scattering channels at vertex 1 and vertex 2, is obtained by putting the two 
respective 2x2 tables on top of one another. For the overall S-channel we thus are 
left with Tab. ffl 




s-ch. 



t-ch. 



u-ch. 



s- ch. 





X 


X 





X 


X 


X 


X 



X 


X 


X 


X 



t-ch. 



X 


X 


X 


X 




X 


X 


X 


X 



X 




X 


X 



X 




X 


X 



X 




X 


X 



X 




X 


X 



u-ch. 



Table 1: Forbidden combinations of energy flow (marked with a X) in all scattering 
channel combinations of vertex 1 and vertex 2 in the overall S-channel. 



Therefore, only six possible combinations remain. Two of them can be excluded 
analytically as we will show in the next section. 
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3.2 Exclusion of configuration Sss 

For the configuration Sss two combinations of energy flow are allowed by Tab. [TJ 
To exclude them, we proceed by a more detailed analysis of the general situation 
expressed by Eqs. (36) and (37) to (40). For our specific cases one has 

2 



a + b 



2ab 



< a b (l — cos ( Zab ) ) < 



a bo{l — cos ( Zab 
1 



< 1 



2a b 



cos 



fZab 



< 



2a b 



The first inequality 



2a b 



< (1 — cos Zab) is always satisfied because 



cos Zab 



and because ao, b > 0. We also have 



u 2 



m 2 
m 2 



(45) 
< 1 

(46) 
(47) 



In addition, four-momentum conservation holds. Thus the a priori eight independent 
entries of the two four-momenta u and v are reduced to four. The on-shellness of 
each mode leads to an additional reduction from four to two: |u| is determined by 
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a, b, and the orientation e u = A = ^ (two angles). One has 



=^ 2 -m 2 =(w-a-5) - m 2 = -2u ( a + 6 ) + 2a6 



- 2 (±\/\u\ 2 + m 2 ^ (a + 6o) + 2ti (a + b) + 2S6 
^±\/ |u| 2 + m 2 ^ (a + &o) = |u| (a + b j e u + a& 



2 



(|u| 2 + m 2 ) (a + 6 

a + b) e u j 2 + (a&y + 2|u| (ab\ (a + b) e u 







u 


2 ( 






|u 


1 



u 



+ ^m 2 (a + 6 j - (ab^j J 
a&j (a + b J e u 



1/2 (s + 6 ) 2 -2((a + b)e u ) 2 



(dt + 6 ) \j-m 2 (a + &o) + + m 2 ((a + b j e u X 

(a + 6 ) - + b) e u ^ 



In the computation of the amplitude only an integration over the orientation e u 
remains when it comes to loop integration. The solutions |u| 1//2 in Eq. (48) must 



be real and positive. This implies that valid configurations satisfy the following 
inequality (argument of square root must be positive): 



-m 2 (a + & ) 2 + (aiy + rh 2 ((a + b) > 
=>- — m 2 (a 2 , + &o + 2a &o) + a 2 ^ 2 , (1 — cos (^ab)) 2 

+ m 2 (a 2 + b 2 + 2a b cos (Zab)) (e u e a+b ) 2 > . (49) 

The values for e u e a+ b, at which the inequality is most likely to be fulfilled, are 
e u e a+b = ±1- For these values the inequality reads as 

afil (1 - cos (Zab)) 2 - 2m 2 h b (1 - cos (Zab)) > . (50) 

To proceed, we need the roots (1 — cos (/ab)) x , 2 of the equation 

a 2 b 2 (1 - cos (Zab)) 2 - 2m 2 h b (1 - cos (Zab)) = . (51) 
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These roots are: 



;i - cos(Zab)) 1 = 



cos (Zab))r 



2m 



, . - 



(52) 
(53) 



Inequality (50) either is solved by values of (1 — cos(Zab)) that are smaller than 



the first solution or larger than the second one because with respect to the variable 
(1 — cos(Zab)) the left hand side of (50) is a parabola with positive curvature. 
(1 — cos (Zab)) 1 < can not be satisfied, and therefore every valid (1 — cos (Zab)) 
has to be larger than 



2ml. 
a b ' 



, , , 2m 2 
(1 - cos(Zab)) > — ^- . 

a b 

We can now compare the two requirements for (1 — cos (Zab)). 



(54) 



From on-shellness and momentum conservation (Eq. (54)): 



i . i \ \ 2m 2 
(1 - cos(Zab)) > . 

a b 



From momentum transfer constraints (Eq. (45)): 

(l_cos(x:ab))<-L- 
2a b 



The upper bound is smaller than the lower bound because fh = 2e > 2y/8n > 1, 
and so we conclude that there are no configurations satisfying momentum transfer 
constraints, on-shellness and momentum conservation simultaneously. Thus a com- 
bination of s and s from the two 4-vertices to the overall S-channel is excluded in 
photon-photon scattering. 



3.3 Energy flow: Overall T- and U-channels 

Similar arguments as used for the overall S-channel can be applied to the overall T- 
and U-channels. The following considerations does not depend on the specific values 
of the external momenta but only on the fact that they are on shell. The overall 
U-channel thus is constrained in the same way as the overall T-channel. (One just 
interchanges the external momenta c and d, see Fig. Q). Let us now put forward 
the arguments for the T-channel. 

The exclusion tables for the constraints at vertex 1 can be obtained by referring to 
Sec.[231 



17 




Figure 8: Overall T-channel. The vertex on the left side is labeled by number 1, the 
one on the right side by number 2. 



s-ch: 1 > I (a — v ) 2 | 



d — u 



— >• 





X 


X 


X 



t-ch: 1 > 



a — d 



{v + uf \ -» 



X 






X 



u-ch: 1 > I (a + u) 2 \ 



~\ 2 

v + d 



— >■ 



X 


X 


X 





The constraints at vertex 2 read: 



s-ch: 1 > 



~\ 2 
v + bj 



(u + cf \ -)> 



X 


X 


X 





t-ch: 1 > I (f — u) 2 | 



b — c 





X 


X 





u-ch: 1 > \(v — c) 2 \ 



b — u 



— > 





X 


X 


X 



In analogy to Sec. 3^ we superimpose all combinations of the tables for the two 
vertices to obtain the result shown in Tab. [2l 
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t-ch. 



u-ch. 



s-ch. 



X 


X 


X 


X 





X 


X 


X 





X 


X 


X 



t-ch. 



X 


X 


X 


X 




X 


X 


X 


X 




X 


X 


X 


X 



u-ch. 



X 


X 


X 





X 


X 


X 





X 


X 


X 


X 



Table 2: Forbidden combinations of energy flow in all scattering-channel combina- 
tions of the overall T-channel (and U-channel). 



For the overall T-channel (or for the overall U-channel) 4 possible combinations 
cannot be excluded. 



3.4 Monte-Carlo analysis of remaining cases 

For 4 out of the 36 configurations in each of the overall S-, T- or U-channels no 
analytical exclusion could be performed. These remaining cases thus are treated 
numerically. To do this, a suitable set of non-redundant variables must be defined. 
As we will see, it is possible to parametrize the overall scattering process by the 
variables (referring to Fig. [6]) ao, bo, the energies of the incoming photons, the 
angle a = Zab G [0, 2ir] between their three-momenta, and the angles C G [— f , f] 
and r] G [0, 27r], that are necessary to describe the spatial orientation of one of the 
outgoing particles. Moreover two angles, 9 G [— f, f] and ip G [0, 2tt], are sufficient 
to describe the kinematic state of the internal particles. 
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Figure 9: A visualization of the involved three-momenta (redundant variables). The 
incoming momenta are blue, the internal ones are red, and the outgoing momenta 
are green. 

In the following we relate the involved four-momenta with respect to these parame- 
ters. First, we exploit four-momentum conservation and that the external photons 
are on shell: 



20 



/ a \ 




/ b 

bo cos (a) 
bo sin (a) 

V o 



/ Co \ 

Co cos (C) COS (fj) 
Co cos (£) sin (77) 
c sin (C) 



d = a -\- b — c 



(55) 



Eq. (|55|) implies 




d z = (a + b-c) = 2ab - 2ac - 2bc 
2 (a b - (a + b ) c - a b cos (a) + c (6 cos (a — rj) + a cos (77)) sin (()) 
a &o(l — cos(a)) 



+ ^0 — b cos (a — 77) sin (^) — a cos (77) sin (() 



independently of the over scattering channel S, T, or U. Momenta of the massive 
modes, however, are dependent on S,T, or U. Both internal particles are on shell. 
For the overall S-channel and for given incoming momenta a, b the absolute value of 
loop three-momentum |u| was related to its orientation e u , which depends on 9 and 



if, in Eq. (48). We have 



u 



ab 



be, 



1/2 f _ v 2 

a + b ) - 2 



± 



be, 



d + bo 



-m 2 ( a + b ) + (ab) + m 2 ( ( a + b) e u 



do + b 



a + b ) e u 



(56) 



The expressions for the overall T- and U-channels can be obtained by interchanging 
momenta: The T-channel relates to the S-channel by exchanging a — > c and b — > 
—a, the U-channel relates to the S-channel by a — > d and b — > —a. Clearly, the 
energy m is determined by on-shellness, and the other internal momentum v is 
completely determined by four-momentum conservation. Whether or not a given 
set of parameter values satisfies all constraints is tested by inserting it into Eqs. 



(55) and (56) before, in turn, the 4- vertex constraints are probed at a given value 



of the dimensionless temperature A. Because A, a , b are not bounded from above 
their to-be-tested values need to be limited in the numerical procedure. We require 
A c = 13.867 < A < 100 [S] and a = ^f, b = ^ < 100. Parameter values are sampled 
randomly in a conditioned way, and parameter sets satisfying the constraints are 
counted (for detailed information, see Appendix [5]). 

3.4.1 Typical hit densities 

For the overall S-channel 87 out of 6.144 x 10 10 tested parameter sets satisfied the 
constraints. In the overall T-channel, we found 1110 valid sets out of 6.144 x 10 10 . 
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Because of the similar structure of the overall U-channel, we expect results very sim- 
ilar to the overall T-channel. These numbers represent the average of all remaining 
scattering combinations in the respective overall channel. Also, the distribution of 



valid parameter values was investigated. For example, Fig. [10] depicts the abun- 
dance of the sum of hits over all non-excluded channel combinations as a function 
of A. 



Hits 
200 



150 



100 



50 



In n 



20 30 40 50 60 70 



Figure 10: This histogram shows the distribution of the sum of hits over all tested 
channel combinations as a function of dimensionless temperature A in the range 
from A c = 13.867 to 100. 



Two things are worth pointing out. First, we see that processes with A > 40 are very 
rare. The highest temperature associated with a valid configuration is A = 63.32, and 
this is an extreme outlier. Second, the abundance is rapidly decaying for A < 18.8. 
In fact, a contribution at a temperature smaller than A = 18.15 never was detected. 
Another interesting distribution is the abundance of energies relative to tem- 



perature T, do and bo, as shown in Fig. 11 As anticipated, there is no obvious 



difference between the distributions of do and bo- The constraints seem to imply an 
upper bound of about 40. Configurations with values of d and b above this bound 
are anyway strongly suppressed by the Bose-Einstein distributions associated with 
external, thermalized photons. 
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Figure 11: Comparison the distribution of the two incoming energies relative to the 
temperature T. Hits are summed over all tested channel combinations. 



3.4.2 Algebraic varieties 

To explore the numerical data further, we zoom into the 6-ip plane about valid 
parameter- value combinations (A, a, (3, 7, do and bo). It was numerically not possi- 
ble to resolve the associated varieties of valid configuration about a particular one 
at once. Therefore we imposed a relaxation of the constraints to broaden the region 
of valid configurations. The relaxation is implemented by virtue of a softening fac- 

(The 



(31) to rt36 



tor T implemented as 1 — > T on the right hand sides in Eqs. 
momentum transfer must only be smaller than T |0| 2 in the relaxed as opposed to 
the physical situation.) For the two shown combinations, Figs. 13 and 14, which are 



representative of all other unexcluded combinations, there are series of four region 
plots. The first (top left) depicts a region of interest of size A9 = 0.4 and Acp = 0.8 
which is centered around the pivotal hit detected during the Monte Carlo test. This 
region is blown up into the second plot (top right) where, in turn, a centered region 
of interest of size AO = 0.04 and A(p = 0.04 is shown. Again, the latter is blown 
up into the third region plot (bottom left). The region of interest marked here and 
shown in full size in plot four (bottom right) has an extent of Ad = 0.004 and 
A(p = 0.004. Values of the softening factor T are chosen to point out the nature of 
the actual, physical variety. The series of region plots in Figs. 13 and 14 have in 
common that the last plot represents this physical situation with T = 1. We need to 
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distinguish the two realizations for each channel combination (corresponding to the 
± in the expression for |u| in Eq. (56)). The solution corresponding to + is labeled 
by the index 1, and the solution corresponding to the minus sign has the index 2. 
The contour plots in the 9-ip plane signal the different constraints by a color code 
as explained in the caption of Fig. |T2 




Figure 12: This region plot gives an color coded representation of the effects of single 
and combined constraints. The varieties marked in orange are associated with real 
and positive solutions |u| to Eq. (56). Varieties, where the momentum transfer 
constraints in the first vertex are satisfied, are marked in green and those, where 
the momentum transfer constraints are fulfilled at the second vertex, are indicated 
in blue. The finally valid variety is represented by the intersection of these three 
varieties. 
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Stti Softening Factor: 150 Stti Softening Factor: 30 




0.15 0.16 0.17 0.18 0.167 0.168 0.169 0.170 



e o 

Figure 13: A visualization of the constraints in the 9-tp plane in the overall S- 
channel with scattering channels t and t for |u| 1 and the following parameter values 
of the pivotal configuration: A =30.9445, a =38.5323, b =47.9434, a =0.985844, 
7] =0.21677, C =1.22281. 
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Figure 14: A visualization of the constraints in the 9-tp plane in the overall T- 
channel with scattering channels u and t for |u| 1 and the following parameter values 
of the pivotal configuration: A =19.2477, a =63.8255, b =0.153577, a =1.46316, 
7] =0.25099, C =0.87816. 



3.5 Bose-Einstein suppression 

Let us now, based on four-momentum conservation, analytically discuss the unex- 
cluded four scattering channel combinations of the overall S-channel. As we have 
seen, the momentum transfer constraints alone were insufficient to exclude certain 
combinations. However, they leave only a small range for possible external energies. 



Because this calculation is generally valid we go back to the nomenclature of Sec. 2.3 
where only massive and massless modes are distinguished: 



V 
R 2 



--0. 



-m 



(57) 
(58) 
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The momentum constraint for cases with a relative minus between the momenta of 
the heavy and the massless modes reads as 



m 2 =f 2p 



R 



+ fh? ± po 



R 



(p-R 

cos ( Z(pR) 



< 1 



< 1 



(59) 



Here =F and the ± represent the energy flow. As one can see from Tab. [TJ the only 
valid combination of energy flow for the four combinations t and t, t and u, u and 
t, or u and u is «o > 0; vo < 0. Comparing this with the constraints of Sec. 3A_ we 
deduce that the following structure is universal: 



m 2 — 2p 



R 



+ m 2 + p 



R 



cos ( Z(pR) 



< 1 . 



(60) 



The absolute-value prescription in (|60|) corresponds to two inequalities as 

cos ( Z(pR)) < 1 . 



—1 < m 2 — 2po 



R 



+ m 2 + po 



R 



(61) 



To eliminate the dependence on Z(pR), we only consider the values =F for the cosine 
for which each respective inequality is most likely satisfied: 



m 2 — 2p 



m 2 — 2p 



R 



_|_ f^2 _ p Q 



R 



+ m 2 + po 



< 1, 
> -1 . 



(62) 
(63) 



It is now straight-forward to deduce the maximally allowed region for po as 

m 2 — 1 fh 2 + 1 

= <Po< 



R 



+ m 2 + 2 



R 



(64) 



R 



+ m 2 - 2 



R 



In the concrete cases of t and t or t and u (64) specializes as 

fh 2 + 1 



m 2 — 1 



2-v/lul +m 2 + 2|u| 



m 2 — 1 



< d < 



<bo< 



u 



fh 2 — 2 lul 



m 2 + 1 



m? + 2 Ivl 



(65) 



2\ v z + m 2 - 2 v 
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For the other two combinations u and t or u and u the dependence on u and v 
switches: 



m 2 — 1 



m 2 + 1 



<a < 



2 A /|v|" + m 2 + 2|v| 
m 2 — 1 



<6n < 



2 A /|vp + m 2 -2|v| 
m 2 + 1 



2\/ u + m 2 + 2 u 



(66) 



2\ u +m 2 -2 u 



The expression resulting from energy conservation at the first vertex is 

a + b =u + vo- 



(67) 



The sum of the energies of the massive particles naturally is high. The idea is 
to search configurations where the upper bounds on the energies of the external 



particles (right-hands sides of (65) and (66)) do not suffice to create two massive 



particles. Regardless of whether one uses (65) or (66), the corresponding inequality 



is 



rh 2 + 1 



m 2 + 1 



+ 



lul + m 2 — 2 lul 



< 



2\l\\\ 2 + m 2 -2|v| 



u| 2 + rh 2 + \/|v| 2 + m 2 



m 2 + 1 



< 



m 2 + l 



2\/|u| 2 + m 2 -2|u| 
m 2 + 1 



~ I 2 i ~ 2 

u + m z 



2\/\v\ 2 + rh 2 - 2 Ivl 



+ \/ I v I 2 + rh 2 . 



(68) 



Defining f(x) = 2 ^^_ 2x 



\/ x 2 + m 2 , this inequality is re-written as: 
<-/(|v|) • 



(69) 



Sufficient but not necessary for (69) to hold is / (x) < 0. Thus, if a configuration 
satisfies 



m 2 + 1 



u| + m 2 < 



(70) 



2\ u +m 2 - 2 u 



and 



m 2 + 1 



v| 2 + m 2 < 



(71) 



2\ v r + w 2 -2 v 
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then we are assured that it does not carry sufficient energy to produce the massive 
loop particles. Coming back to / (x) < 0, we have 



ra 2 + 1 



V x 2 + m 2 < 



2\/x 2 + fa 2 — 2x 
=>- 2 (x 2 + m 2 ) — 2x \/ x 2 + fa 2 > fh 2 + 1 

> 1 + 2x\/x 2 + fh 2 - 2x 2 - fh 2 < . 



(72) 



Looking at the corresponding equation^ we obtain a bound on the allowed region 
for x 



1 + 2xV x 2 + m 2 — 2x 2 — m 2 = . 



(73) 



This is solved by X1/2 = ±| (fh 2 — 1). Each of these solutions is twofold degenerate. 
Because we are looking at the absolute values of three momenta it follows that 
x > 0. Because fh > 1 this means that only the root with + is admissible. But at 



x 



Um 2 



1) the derivative of function g with respect to x is positive: 



dg_ 

dx 



dg (x) 
dx 



2y/ x 2 + m 2 + 



2a; 2 



Vx^+ 



m 2 



- Ax 



+ 



(m 2 — 1) + m 2 



I {fa 2 - if 
| (fh 2 — l) 2 + fh 2 



2 (fh 2 - 1) > 



Inequality (74) holds because of the following: 



(74) 



1 



I (™2 _ i 

2W-(m 2 -l) 2 + m 2 + 2 

4' /l/^9 t\2 



2 (m 2 — l) z + m 2 



- 2 (fh 2 - l) > 



2 Q (fa 2 - l) 2 + m 2 ^ + ^ (fh 2 - l) 2 > 2 (fa 2 - l) 
<^ (fh 2 - l) 2 + 2m 2 > 2 (m 2 - l) W - (m 2 - l) 2 + fh 2 



[ra 2 — 1) + ra 2 



& (fa 2 - l) 4 + 4m 4 + 4m 2 (fh 2 - if > A (fh 2 - if Q (m 2 - l) 2 + m 2 ^ 



<^ 4m 2 > . 



(75) 



Let us now summarize our result: For non- negative x there is only one location 
x = \ (fh 2 — 1) where g (x) = 0. At this point the slope of function g is positive, see 



'This equation is obtained by setting the left-hand side of (172]) equal to zero 
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(74). Thus g (x) is smaller than zero (as demanded by (72)) for values of x smaller 
than | (m 2 — 1). This implies that for < x < | (m 2 — 1) configurations of external 
photons energies cannot create two massive particles. In other words: |u| and |v| 
have to be at least of the value \ (m 2 — 1) in the energy-conservation equation at 



vertex 1, Eq. (67). That is 



u ,v > 



m 2 — 1) + m 



vi 2 



(76) 



The propagator of a massive mode comes with a factor of % ' 



I ftp I 

e t 



the 



Bose-Einstein distribution. For low temperatures this distribution suppresses the 



process strongly, see Fig. 15 Because of (76) the following inequalities are valid: 
1 1 



exp (u ) - 1 ' exp (v ) 



< 



1 



exp 



2tt 
AV2 



if 



m? 



(77) 



1 



Because m = 2e (e > \/8ir) [20] is a function of the dimensionless temperature A, 
which is numerically available [8J, we may infer the upper bound for the product of 



the two Bose-Einstein distributions as a function of A, see Fig. 15 Obviously, there 




Figure 15: The upper limit for the product of Bose-Einstein distributions as a 
function of dimensionless temperature A (blue area is not excluded). The left panel 
shows this function in the domain A c < A < 100 while the right plot depicts it in 
the domain X r < A < 30. 



is strong suppression near the critical temperature due to fa — > oo for A \ A c . For 
a numerical impression, see also Tab. [3| For A > 40 we only obtain suppression 
factors > 10" 4 . But, as suggested by Monte-Carlo sampling, for temperatures this 



high the 4- vertex constraints practically are never satisfied, see Sec. |3.4.1j On the 
other hand, for A < 40 the hit densities are sizable but the Bose-Einstein suppression 
is severe. 



30 



A 


14 


15 


20 


max (n B (uo) x n B {v )) 


3.7* 10~ 159 


2.3 * 10" 60 


1.4* 10~ 18 


A 


25 


30 


40 


max (n B (uo) x «b(^o)) 


5.6 * KT 11 


1.0 * 10" 7 


8.6 * 10" 5 


A 


50 


100 


1000 


max(n B (?io) x ^(vq)) 


1.8* 10~ 3 


3.0 * 10" 1 


9.8* 10 2 



Table 3: Examples for upper bounds on the product of the two Bose-Einstein distri- 
butions, associated with the massive modes in the loop, which occur in the integrand 
of the amplitude for photon-photon scattering. 



4 Summary and Conclusions 

In the present work we have investigated systematically how the unitary-gauge con- 
straints of the effective theory for deconfming SU(2) Yang-Mills thermodynamics 
severely limit the contributions to the amplitude for one-loop photon-photon scat- 
tering. Only one type of Feynman diagrams with two 4-vertices is admissible to 
mediate this process, and a large part of channel and energy-sign combinations for 
the scattering through these vertices is analytically excluded relying on a partial 
set of all energy-flow and momentum-transfer constraints. However, out of 108 
scattering-channel and energy-sign combinations 12 cannot be excluded in this way. 
These remaining cases were analyzed by Monte Carlo sampling subject to all con- 
straints generating hit densities which decay very rapidly with temperature. We 
have also investigated the admissible algebraic varieties in the vicinity of selected, 
pivotal Monte-Carlo hits to demonstrate how filamented these regions are. On the 
other hand, we have shown by analytical arguments for the overall S channel that for 
small temperatures the product of Bose-Einstein distributions, associated with the 
thermal on-shell propagation of the massive loop modes, suppresses the amplitude 
very strongly at small temperatures. An analogous analytic access to the overall T 
and U channels could not be exploited but we expect similar behavior also in these 
channels. All in all, our results suggest that photon-photon scattering within the de- 
confining SU(2) Yang-Mills plasma practically does not occur. The optical theorem 
thus would exclude the annihilation of two massive modes into two photons or the 
pair creation of massive modes by photons which is in agreement with experiment^ 
A higher dimensional generalization (more than two 4-vertices per diagram) of 
the technology of energy-flow exclusion developed in the present work may be key 
to proving the termination of the expansion of thermodynamical quantities into 

10 In [19] it is explained at length why the gauge-group factor U(l)y- of the present Standard 
Model of Particle Physics should be interpreted as the Cartan group of dynamically broken SU(2) 
Yang-Mills theory of scale A ~ 10~ 4 eV. 
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irreducible bubble diagrams at a finite loop order, conjectured in [27] and based on a 
counting of constraints versus independent loop-integration variables in dependence 
on loop order. 



5 Appendix 



Here we sketch the sampling strategy referred to in Sec. |3.4| Because we would 
like to do justice to the different scales of the sampling ranges and since runtime 
resources were limited we nested the sampling of random variables and required 
that noncompact variables are tested much more often than compact ones. The 
algorithm is sketched below in terms of a pseudocode: 



repeat 60 times: 

A = random real number G [13.867, 100]; 
repeat 4 times: 

a = random real G [0, 2ir]; 
repeat 4 times: 

t] = random real G [0, 2tt]; 
repeat 2 times: 

( = random real G [— f, fl; 
repeat 2 times: 

9 = random real G [— | , f] ; 
repeat 4 times: 

ip = random real G [0, 2ir]; 
repeat 10 4 times: 

a = random real G [0, 100]; 
b = random real G [0, 100]; 
test if constraints for given channel-combination 
are satisfied for the given parameters 
A, a , S , a, T], C, 0, tp; 
if yes: 

increase hitcount by one; 
save parameters; 
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